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The Cooling Storage Ring external-target experiment (CEE) spectrometer is used to study the nuclear matter 
created in heavy-ion collisions at \/Syx = 2.1-2.4 GeV with the aim to reveal the quantum chromodynamics 
phase structure in the high-baryon density region. Collective flow is considered an effective probe for evaluating 
the properties of media during high-energy nuclear collisions. One of the main functions of the zero-degree 
calorimeter (ZDC), a subdetector system in the CEE, is to determine the reaction plane in heavy-ion collisions. 
This step is crucial for measuring the collective flow and other reaction-plane-related analyses. In this paper, 
we illustrate the procedures for event-plane determination using the ZDC. Finally, isospin-dependent quantum 
molecular dynamics model-based predictions of the rapidity dependence of the directed and elliptical flows for 
p, d, t, 3He, and “He, produced in 2.1 GeV U+U collisions, are presented. 
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I. INTRODUCTION 


Quantum chromodynamics (QCD) predicts a phase tran- 
sition from hadronic to deconfined partonic matter at suffi- 
ciently high temperatures and/or densities [1]. Results from 
top RHIC and LHC energies indicate that a new form of mat- 
ter with a low viscosity and high temperature, quark-gluon 
plasma (QGP), has been produced [2-7]. Lattice QCD cal- 
culations predict that the phase transition from hadronic mat- 
ter to the QGP phase is a smooth crossover in the vanishing 
baryon chemical potential (ug) region [8]. A first-order phase 
transition is expected in a finite baryon chemical potential re- 
gion, and thus, determining the phase structure of the QCD is 
a major research goal in the field of medium- and high-energy 
heavy-ion collisions [9-12]. 

The Cooling Storage Ring external-target experiment 
(CEE) is a spectrometer that is employed to investigate the 
properties of nuclear matter in the 2.1-2.4 GeV energy region 
in the center-of-mass frame [13]. The CEE primarily allows 
near-full-space measurements of charged particles in heavy- 
ion collisions and provides experimental data for studying im- 
portant scientific problems, such as spin- and isospin-related 
nuclear forces, nuclear matter equations of state, and QCD 
phase structures at high baryon number densities [14-16]. 
This offers valuable research opportunities for QCD phase 
diagram studies in low-temperature and high-baryon-density 
regions. 

The event anisotropy of final-state particles relative to the 
reaction plane in momentum space, also known as collective 
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flow [17], is important for evaluating the properties of me- 
dia created in heavy-ion collisions. Flow coefficients, such 
as directed flow vı and elliptical flow v2, are characterized 
by harmonic coefficients in the Fourier expansion of the az- 
imuthal distribution of the final particles with respect to the 
reaction plane. In heavy-ion collisions, the driving force 
of the collective flow is the initial anisotropy in coordinate 
space. It rapidly diminishes as a function of time, and this 
phenomenon is known as the self-quenching effect. Thus, the 
collective flow is sensitive to the details of the expansion of 
nuclear matter during the early collision stage. The directed 
flow vı was predicted to be sensitive to the effective equation 
of state [18-20]. The elliptic flow v2 is sensitive to the con- 
stituent interactions and degrees of freedom [21-23]. CEEs 
can reveal the collective flow in heavy-ion collisions at \/5.y 
= 2.1-2.4 GeV. This will help us study the medium prop- 
erties and further search for possible QCD phase transition 
signals [24-26]. One of the principle functions of the zero- 
degree calorimeter (ZDC), a subdetector of the CEE, is to de- 
termine the reaction plane in nuclear—nuclear collisions. The 
reconstructed reaction plane (usually called the event plane) is 
crucial for many measurements, such as collective flow [27— 
29], azimuthal HBT [30], and CME-related observables [31— 
34]. 

In this paper, we introduce the necessary acceptance cor- 
rections and calibrations for event-plane determination from 
the CEE-ZDC. Furthermore, the isospin-dependent quantum 
molecular dynamics (IQMD) model [35] is used to predict 
collective flow from a typical CEE energy (Syn = 2.1 GeV). 


Il. CEE-ZDC 


Figure 1(a) shows a schematic of the CEE spectrometer. 
The detector subsystem consists of a superconducting dipole 
magnet used to deflect charged particles; a silicon pixel posi- 
tioning detector (SiPiX, Beam Monitor) to measure the posi- 
tion and time of the incident beam as well as the primary col- 
lision vertex [36]; A time projection chamber (TPC) is used 


to reconstruct particle trajectories and identify particles [37], 
and a time-of-flight chamber (TOF) is employed to extend 
particle identification to high momentums (p > 2 GeV/c) 
The TOF chamber contains a start-time detector (TO) [38], 
an inner time-of-flight detector (TOF) [39], and an end-cap 
time-of-flight detector (eTOF) [40]. In addition, multi-wire 
drift chambers (MWDCs) are designed to track charged par- 
ticles at forward rapidity and identify particles via momen- 
tum measurements [41]. The ZDC is used to measure the 
patterns (deposited energy and incident position) of forward- 
going charged particles emitted from nuclear—nuclear colli- 
sions [42]. 

The ZDC is proposed to be installed behind all the other 
subdetectors. The beam direction is defined as the positive 
Zaxis, and the ZDC is located at Z = 295-299 cm, facing 
the original incident beam direction; ilts geometry is shown 
in Fig. 1(b). The ZDC detector cross-plane is a wheel with 
radius R ranging from 5 to 100 cm, and the vacuum pipe 
carrying the nuclear beam passes through the inner hole of 
this ZDC wheel. It consists of 24 sectors that subtend an az- 
imuthal angle of 15°. Each sector is divided into eight mod- 
ules forming eight rings in the full ZDC plane. The sensi- 
tive volume of the ZDC is composed of a plastic scintillator, 
and the current design uses the BC-408 material from Saint- 
Gobain [43]. The photons are produced inside the scintil- 
lator through the deposited energy of the incident particles 
and are then transported through a plastic light guide into the 
quartz window of a traditional PMT. ZDC covers the pseudo- 
rapidity range between 1.8 and 4.8, allowing the determina- 
tion of the centrality and event plane in the forward rapidity 
region and minimizing autocorrelations from middle rapidity 
analyses [17, 44]. 
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Fig. 1. (Color online) (a) Schematic of the CEE detector. (b) 
Schematic of the ZDC detector. 


I. EVENT PLANE DETERMINATION FROM THE 
CEE-ZDC 


In the study of the event plane, the simulation input of 735U 
+ 238 collisions at 500 MeV/u was obtained from the IQMD 
generator [35]. The IQMD model was developed based on the 
quantum molecular dynamics (QMD) model [45] consider- 
ing isospin effects. The detector environment was simulated 
using GEANT4 [46]. One million IQMD simulated events 


were generated in the range of the nuclear impact parameter 
b, which is the transverse distance of the projectile from the 
target nucleus, 0 < b < 10 fm, with 0.1 million events for 
each b interval of 1 fm. 

The reaction plane in the nucleus-nucleus collision is de- 
fined by the vector of the impact parameter and beam direc- 
tion. As the impact parameter could not be directly measured 
in the experiment, the reaction plane was estimated using the 
standard event plane method [17, 47]. The first-order har- 
monic event plane W is calculated using the event flow vec- 


tor Qi: 
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(1) 
where the sum exceeds all particles used in the flow vector 
calculation. Quantities ġ; are the azimuths in the laboratory 
frame. The weight w; is defined by the deposited energy AE 
of particle ¿ collected by ZDC detector. As it is related to the 
mass and transverse momentum pr of the particle, whereas 
the pr weight is commonly applied in flow analysis to opti- 
mize the event plane resolution [47]. The smearing effect of 
the deposited energy is considered by Equ. 2 
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where L is the distance from the hit position to the geomet- 
ric center of the sector. h is the charge of the final particle. 
The term 1 — $(+4)? is used to describe the resolution of 
the deposited energy at the edge of the sector. and the term 
8 + 3(h — 8) is used to simulate the saturation effect of the 
deposited energy resolution for heavy nuclei (h > 8) [48]. 

Because the finite multiplicity limits the estimation of the 
reaction plane, it yields a resolution factor R which is defined 
by Equ. 3. In this study, we focus on the first-order harmonic 
event plane because vı is more significant than higher-order 
flows in the collision energy range covered by CEE. 


R, = (cos (Wi ep — Vi RP)) (3) 


The magnetic field direction was perpendicular to the beam 
direction at CEE. Thus, the charged particles in the final state 
are deflected by the magnetic field and hit one side of the ZDC 
detector more, as shown in Fig. 2(a) Owing to the asym- 
metric ZDC acceptance, The reconstructed event-plane angle 
was not isotropic within the laboratory frame. but is biased 
towards the 7-azimuth. The acceptance bias caused by the 
magnetic field introduces additional nonphysical anisotropy 
for the detected collision events, and this effect should be re- 
moved as it distorts the event plane reconstruction. Therefore, 
we introduce a position weight to calibrate asymmetric accep- 
tance. 

The core idea of the position weight is to correct the asym- 
metric acceptance of the ZDC caused by the magnetic field. 
Owing to the deflection of charged particles in the magnetic 
field, The left side of the ZDC detector receives more hits. We 


Fig. 2. (Color online) (a) Hit distribution obtained from the ZDC 
with a collision impact parameter of 5 < b < 6 fm. (b) Hit distribu- 
tion obtained from the ZDC after position weight correction with a 
collision impact parameter of 5 < b < 6 fm. 


assigned a weight P of less than | to the hits on the left side 
to correct this effect. The weight was calculated based on a 
two-dimensional X — Y hit distribution, as defined in Equ. 
4 is the ratio of the number of hits on the right side to those 
on the left side. In addition, the deposited energy AE was 
used as the weight when calculating the number of hits. be- 
cause it is related to particle mass. It can be observed that the 
acceptance of ZDC is symmetric after applying the position 
weight, as shown in Fig. 2(b). 
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Fig. 3. (Color online) Event plane distributions before position 
weight correction (black line), after position weight correction (red 
line), and after position weight + shift corrections (blue line). 


The black line in Fig. 3 shows the event plane distribution 
before position weight correction. With an ideal detector, the 
event plane distribution should be flat because the possible 
direction of the impact parameter b is random in the 27 az- 
imuths of the transverse plane in the laboratory frame. It is not 
flat but peaks around Y, ~ 7 owing to the asymmetric accep- 
tance of ZDC as discussed above. Correspondingly, one can 


see that the resolution difference between the left (azimuth 
of the reaction plane: 1/2 to 37/2) and right sides (—7/2 to 
1/2) of ZDC is significant, as shown in Fig. 4(a). After ap- 
plying the position weights defined in Equ. 4: The unflatness 
of the event plane is significantly reduced, as indicated by the 
red line in Fig. 3. The resolution difference between the left 
and right sides of ZDC is significantly reduced, as shown in 
Fig. 4(b). This indicates that the position weight naturally 
corrected the acceptance asymmetry of ZDC. 
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Fig. 4. (Color online) (a) Resolution of 1% order event plane as a 
function of impact parameter b without position weight. (b) Resolu- 
tion of 1% order event plane as a function of b with position weight. 
(c) Resolution of 1% order event plane as a function of b with posi- 
tion weight and shift correction. 


The event plane distribution is not perfectly flat after the 
position weight, as shown in Fig. 3. Consequently, the reso- 
lution difference between the left and right sides of ZDC was 
still visible. Therefore, the shift method is used to force the 
event plane to be flat [47]. A shift angle AW, is applied to 
correct the event plane. and AW, is calculated event-by-event 
using the following equation: 


Y =U, + AV, 
20 
AY, = 5 z [— (sin (iW,)) COS (iW) + (cos (iW,)) sin (iv,)] 


i=l 
(5) 
where the brackets refer to the average over the events in 
the same centrality bins. W, is the position-weight-corrected 
event plane azimuth, and Ų/ is the event plane angle with 
shift calibration. After the shift calibration, a flat event plane 
distribution was achieved, as indicated by the blue line in Fig. 
3. The resolution between the left and right sides was consis- 
tent, as shown in Fig. 4(c)). 

In the experiment, the event plane calculated from differ- 
ent rapidity windows helped us understand the systematic un- 
certainties in the flow measurements. Correspondingly, the 
event planes from ZDC sub-rings, which correspond to dif- 
ferent rapidity windows, were studied. Figure 5 shows the 
1% order event plane resolution from ZDC sub-ring radius 
52.5 < R < 76.25 cm without a position weight 5(a), with a 
position weight 5(, ) and with a position weight and shift cor- 
rection 5(c). These results indicate that the position weight 
and shift methods work well for the event plane calculated by 
ZDC subring. 

After eliminating the resolution difference due to asym- 
metric acceptance using the position weight and shift meth- 
ods, the 1**-order event plane resolution from the ZDC was 
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Fig. 5. (Color online) (a) Resolution of 1% order event plane, de- 
termined from the ZDC sub-ring (radius: 52.5 < R < 76.25 cm) 
without position weight. (b) Resolution of 1% order event plane, de- 
termined from the ZDC sub-ring (radius: 52.5 < R < 76.25 cm) 
with position weight. (c) Resolution of 1% order event plane, deter- 
mined from the ZDC sub-ring (radius: 52.5 < R < 76.25 cm) with 
position weight and shift corrections. 


calculated using the two-sub-event plane method [17]. The 
full event was divided randomly into two independent sub- 
events with equal tracks, and the event-plane resolution was 
estimated by correlating the two sub-events. as defined by 
Eq. 6: 


Ri sa» = y (cos (Uf — WP) ) 
= VT/2xe 0’? (Io (xX?/2) + L (x?/2)) 


where A and B denote the two subevents. As x is proportional 
to the square root of the multiplicity and a full event with 
two particles as subevents, The full event-plane resolution is 
obtained as follows: 


(6) 


Ran = R (V2xsw ) o 


The resolution of 15%t-order event plane as a function of im- 
pact parameter, determined from the ZDC whole ring, is com- 
pared with the 1**-order event plane resolution, determined 
from the STAR event plane detector, for Au+Au collisions at 
J/Sxn = 3.0 GeV [28] in Fig. 6. The event plane resolution 
from CEE-ZDC reached ~ 90% in the middle-central colli- 
sions (4 < b < 7 fm). The resolution of ZDC is better in the 
region of b < 4 fm but worse for b > 4 fm, which is proba- 
bly due to the different sizes of the gold and uranium nuclei, 
experimental acceptance, and detector performance. 

We also systematically investigated the effects of the ZDC 
detector thickness, hit efficiency, energy resolution, and 
model dependence on the first-order event plane resolution. 
As shown in Fig. 7, where the solid red dots represent the de- 
fault conditions: a ZDC thickness of 4 cm, hit efficiency of 
100%, and default energy smearing, as shown in Eq. 2 and 
heavy nuclei from IQMD generator de-excitation. The ef- 
fects of the different variables were investigated individually. 
The resolution of 1%t-order event plane decreases slightly 
as the ZDC detector thickness decreases, as shown in Fig. 
7(a). This is because a more accurate measurement of the 
deposited energy is achieved with a thicker ZDC. Fig. 7(b) 
shows the hit efficiency dependence of 1*t-order event-plane 
resolution. The ZDC hit efficiency was reduced to 90%, and 
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Fig. 6. (Color online) Resolution of 1% order event plane, deter- 
mined from the ZDC whole ring, as a function of impact parameter 
compared with that of 1% order event plane determined from STAR 
event plane detector for Au+Au collisions at \/S5yx = 3.0 GeV. 
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Fig. 7. (Color online) (a) Effect of ZDC thickness on the 1% order 
event plane resolution. (b)Effect of ZDC hit efficiency on the 1** 
order event plane resolution. (c)Effect of ZDC energy resolution on 
the 1% order event plane resolution. (d)Effect of heavy nuclei de- 
excitation on the 1% order event plane resolution. 


the event-plane resolution remained almost unchanged. The 
effect of ZDC energy resolution is investigated by applying 
additional Gaussian smearing to the deposited energy, where 
Gaussian(1, 0.5) has a center value of 1 and width of 0.5, 
and Gaussian(1, 1) has a center value of 1 and width of 1. A 
smaller Gaussian width indicates a better energy resolution. 
As the energy resolution decreases, the first-order event-plane 
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Fig. 8. (Color online) vı and v2 as functions of rapidity, for protons, deuterons, tritons, 3He, and *He with a collision impact parameter 
1 < b < 4 fm, derived from IQMD 7°8U + ?3SU 500 MeV/u (, /Sxx = 2-1 GeV). The v1 slopes are extracted using the equation: y = az+bz?. 
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Fig. 9. (Color online) Atomic number A scaled vı slope (upper 
panel) and v2 (lower panel) at middle rapidity for protons, deuterons, 
tritons, 3He, and “He determined from HADES, STAR experiments, 
and IQMD model calculations for CEE. The brackets represent the 
systematic uncertainties in the experimental data. 


resolution decreases by approximately 5-10%, as shown in 
Fig. 7(c). Figure 7(d) shows the relationship between the 
ZDC event plane resolution and the IQMD heavy nuclei de- 
excitation, where ”out”/”in” means the heavy nuclei are de- 
excitation or not. The resolution estimated using the IQMD 
model with heavy nuclei de-excitation was slightly higher 
than IQMD without heavy nuclei de-excitation because the 
multiplicity was higher in the former case. 


IV. COLLECTIVITY FLOW PREDICTIONS FROM IQMD 
MODEL 


Collective flow is sensitive to the details of the expansion of 
the medium produced during the early collision stage. Flow 
measurements at CEE provide information on the QCD phase 
structure in the high-baryon-density region. Collectivity flow 
predictions were presented for a typical CEE based on IQMD 
model. Figure 8 shows vı and v2 as functions of rapidity for 
protons, deutons, tritons, *He, and He with an impact param- 
eter of 1 < b < 4 fm from IQMD 72°U + ?38U collisions at 
500 MeV/u (Syn = 2.1 GeV). The vı slope values extracted 
using y = ax + bz? strongly depend on the number of nu- 
clei. The və values are negative in the middle rapidity owing 
to the squeeze-out effect; the medium expansion is shadowed 
by spectator nucleons, and particles are preferentially emit- 
ted in the direction perpendicular to the reaction plane [17], 
whereas v2 becomes positive in the forward rapidity as the 
squeeze-out effect becomes weak. Similar to vı slope, v2 
values showed a strong dependence on the number of nuclei. 

Figure 9 presents the d(v;/A)/dy and v2 for pro- 
tons, deuterons, tritons, *He, and He determined from 
HADES [49] and STAR [50] experiments together with the 
IQMD model calculations 3. where A denotes the atomic 
number. v;/A represents the directed flow carried by each 
nucleon in the light nuclei, and the scaling behavior suggests 
a coalescence production mechanism of the light nuclei dur- 
ing heavy-ion collisions. və is calculated in the rapidity range 
of —0.1 < y < 0 for the STAR experiment and IQMD model 
calculations and —0.05 < y < 0.05 for the HADES experi- 
ment. The atomic-number-scaled vı slope from HADES and 


3 Unlike in the experiment, the centrality here is determined from the impact 
parameter in the model calculations. 


IQMD showed a decreasing trend with an increase in atomic 
number, whereas the STAR data weakly depend on the atomic 
number in the collisions at ,/s,,, = 3 GeV. The absolute value 
of v2 from HADES decreased with increasing atomic num- 
ber, whereas the results for STAR and IQMD remained al- 
most unchanged with atomic number. This may indicate that 
light nuclei are not purely formed by the coalescence mech- 
anism in Au+Au collisions at \/5,,, = 2.4 GeV, Coalescence 
was the dominant production mechanism for Au+Au ,/5,, = 
3.0 GeV. The production of light nuclei in the IQMD model 
is a mixture of light nuclei fragments and the coalescence of 
nucleons and light nuclei. The dominant production mecha- 
nism in the IQMD model depends on the collision energy and 
parameter settings. The predictions provided by the IQMD 
model in U + U collisions at \/s,,, = 2.1 GeV will be vali- 
dated in future CEEs. 

Future measurements of vı and v2 will help us to examine 
the equation of state of the produced nuclear matter at CEE 
energies [51, 52] and understand the production mechanism 
of light nuclei in the high baryon density region [15, 53—56]. 


V. SUMMARY 


In this paper, we elucidate the procedures for event plane 
determination from the ZDC at the CEE. The calculated val- 
ues determined using the IQMD Monte Carlo event generator 
(500 MeV/u 735U + ?38U) were incorporated as inputs, and 
the detector environment was simulated using GEANT4. 

To correct for the bias caused by the dipole magnet, a 
position-dependent weight was introduced to calibrate the 
asymmetric acceptance. After an additional shift correction, 
an outstanding first-order event plane resolution of ~ 90% 
was obtained for middle-central collisions (4 < b < 7 fm). 
Herein, the collective flows vı and v2, as functions of rapidity 
for p, d, t, 3He, and He in middle central collisions, are pre- 


[1] P. B. Munzinger and J. Stachel, The quest for the quark- 
gluon plasma. Nature 448, 302-309 (2007). doi: 10.1038/na- 
ture06080 

[2] I. Arsene et al., Quark gluon plasma and color glass 
condensate at RHIC? Perspective of the BRAHMS 
experiment. Nucl. Phys. A 757, 1-27 (2005). doi: 
10.1016/j.nuclphysa.2005.02.130 

[3] B. B. Back et al., The PHOBOS perspective on discov- 
eries at RHIC. Nucl. Phys. A 757, 28-101 (2005). doi: 
10.1016/j.nuclphysa.2005.03.084 

[4] K. Adcox et al., Formation of dense partonic matter in rela- 
tivistic nucleus-nucleus collisions at RHIC: Experimental eval- 
uation by the PHENIX collaboration. Nucl. Phys. A 757, 184- 
283 (2005). doi: 10.1016/j.nuclphysa.2005.03.086 

[5] J. Adams ef al., Experimental and theoretical challenges 
in the search for the quark gluon plasma: The STAR 
Collaboration’s critical assessment of the evidence from 
RHIC collisions. Nucl. Phys. A 757, 102-183 (2005). doi: 
10.1016/j.nuclphysa.2005.03.085 


sented based on the IQMD model. These results were com- 
pared with the experimental data obtained from 2.4 GeV and 
3 GeV Au+Au collisions in the HADES and STAR experi- 
ments, respectively. 

The measurements from the HADES and STAR experi- 
ments suggest that coalescence is the dominant production 
mechanism of light nuclei at 3 GeV, whereas light nuclei frag- 
mentation and coalescence are both important at 2.4 GeV. The 
predictions of the IQMD model at 2.1 GeV will be validated 
in future CEEs. 


ACKNOWLEDGMENTS 


We thank Prof. Li Ou and Zhigang Xiao for generating the 
IQMD data and for fruitful discussions. 


AUTHOR CONTRIBUTIONS 


All authors contributed to the study conception and de- 
sign, material preparation, data collection and analysis. The 
first draft of the manuscript was written by Li-Ke Liu and 
Shu-Su Shi, and all authors commented on previous versions 
of the manuscript. All authors read and approved the final 
manuscript. 


DATA AVAILABILITY STATEMENT 


The data that support the findings of this 
study are openly available in Science Data Bank 
at https://doi.org/10.57760/sciencedb.09 132 and 
https://cstr.cn/3 1253.1 1.sciencedb.09132. 


[6] A. Bazavov et al., The chiral and deconfinement aspects of 
the QCD transition. Phys. Rev. D 85, 054503 (2012). doi: 
10.1103/PhysRevD.85.054503 

[7] K. Fukushima and C. Sasaki, The phase diagram of nuclear and 
quark matter at high baryon density. Prog. Part. Nucl. Phys 72, 
99-154 (2013). doi: 10.1016/j.ppnp.2013.05.003 

[8] Y. Aoki, G. Endrodi, Z. Fodor et al., The Order of the quantum 
chromodynamics transition predicted by the standard model of 
particle physics. Nature 443, 675-678 (2006). doi: 10.1038/na- 
ture05120 

[9] K. Fukushima and T. Hatsuda, The phase diagram of 
dense QCD. Rept. Prog. Phys. 74, 014001 (2011). doi: 
10.1088/0034-4885/74/1/014001 

[10] A. Bzdak et al., Mapping the phases of quantum chromody- 
namics with beam energy scan. Phys. Rept. 853, 1-87 (2020). 
doi: 10.1016/j.physrep.2020.01.005 

[11] X. Luo, S. Shi, N. Xu et al., A study of the properties of the 
QCD phase diagram in high-energy nuclear collisions. Particle 
3(2), 278-307 (2020). doi: 10.3390/particles3020022 


[12] H. Huang et al., Collective excitation in high-energy nuclear 
collisions — In memory of professor Lianshou Liu. Symmetry 
15 499 (2023). doi: 10.3390/sym15020499 

[13] L. Lü et al., Conceptual design of the HIRFL-CSR external- 
target experiment. Sci. China Phys. Mech. Astron. 60(1), 
012021 (2017). doi: 10.1007/s11433-016-0342-x 

[14] C. J. Horowitz et al., A way forward in the study of the sym- 
metry energy: experiment, theory, and observation. J. Phys. G 
41, 093001 (2014). doi: 10.1088/0954-3899/41/9/093001 

[15] S. Zhang et al., Searching for onset of deconfinement via hy- 
pernuclei and baryon-strangeness correlations. Phys. Lett. B 
684, 224-227 (2010). doi: 10.1016/j.physletb.2010.01.034 

[16] A. Andronic et al., Hadron production in ultra-relativistic nu- 
clear collisions: Quarkyonic matter and a triple point in the 
phase diagram of QCD. Nucl. Phys. A 837, 65-86 (2010). doi: 
10.1016/j.nuclphysa.2010.02.005 

[17] S. A. Voloshin, A. M. Poskanzer, R. Snellings, Collective phe- 
nomena in non-central nuclear collisions. Landolt-Bornstein. 
23, 293-333 (2010). doi: 10.1007/978-3-642-01539-710 

[18] S. A. Bass et al., Microscopic models for ultrarelativistic heavy 
ion collisions. Prog. Part. Nucl. Phys. 41, 255-369 (1998). doi: 
10.1016/S0146-6410(98)00058-1 

[19] J. Steinheimer et al., The high-density equation of state in 
heavy-ion collisions: constraints from proton flow. Eur. Phys. J. 
C, 82(10), 911 (2022). doi: 10.1140/epjc/s 10052-022-10894-w 

[20] D. Oliinychenko, A. Sorensen, V. Koch et al., Sensitivity of 
Au+Au collisions to the symmetric nuclear matter equation of 
state at 2 — 5 nuclear saturation densities. arXiv:2208.11996 

[21] L. Adamczyk et al., Centrality and transverse momentum de- 
pendence of elliptic flow of multistrange hadrons and ¢ me- 
son in Au+Au collisions at 4/Syy = 200 GeV. Phys. Rev. Lett. 
116(6), 062301 (2016). doi: 10.1103/PhysRevLett.116.062301 

[22] L. Adamczyk et al., Measurement of D? azimuthal anisotropy 
at midrapidity in Aut+Au collisions at \/Syxy = 200 GeV. 
Phys. Rev. Lett. 118(21), 212301 (2017). doi: 10.1103/Phys- 
RevLett.118.212301 

[23] S. Shi, An experimental review on elliptic flow of strange 
and multistrange hadrons in relativistic heavy ion colli- 
sions. Adv. High Energy Phys. 2016, 1987432 (2016). doi: 
10.1155/2016/1987432 

[24] Y. Nara and A. Ohnishi, Mean-field update in the JAM micro- 
scopic transport model: Mean-field effects on collective flow 
in high-energy heavy-ion collisions at ,/Syy = 2-20 GeV ener- 
gies. Phys. Rev. C 105(1), 014911 (2022). doi: 10.1103/Phys- 
RevC.105.014911 

[25] Y. Nara, A. Jinno, K. Murase et al., Directed flow of A 
in high-energy heavy-ion collisions and A potential in dense 
nuclear matter. Phys. Rev. C 106(4), 044902 (2022). doi: 
10.1103/PhysRevC.106.044902 

[26] S. Lan and S. Shi, Anisotropic flow in high baryon density re- 
gion. Nucl. Sci. Tech. 33(3), 21 (2022). doi: 10.1007/s41365- 
022-01006-0 

[27] L. Adamczyk et al., Beam-energy dependence of the directed 
flow of protons, antiprotons, and pions in Au+Au collisions. 
Phys. Rev. Lett. 112(16) 162301 (2014). doi: 10.1103/Phys- 
RevLett.112.162301 

[28] M. S. Abdallah et al., Disappearance of partonic collectivity in 
J/Sxw = 3 GeV AutAu collisions at RHIC. Phys. Lett. B 827, 
137003 (2022). doi: 10.1016/j.physletb.2022.137003 

[29] Y. Nara et al., Dynamical evolution of heavy-ion collisions. 
doi: 10.1007/978-981-19-4441-33 

[30] L. Adamczyk et al., Beam-energy-dependent two-pion inter- 
ferometry and the freeze-out eccentricity of pions measured in 


heavy ion collisions at the star detector. Phys. Rev. C 92 014904 
(2015). doi: 10.1103/PhysRevC.92.014904 

[31] Fukushima et al., The chiral magnetic effect. Phys. Rev. D 78, 
074033 (2008). doi: 10.1103/PhysRevD.78.074033 

[32] L. Adamczyk et al., Search for the chiral magnetic effect via 
charge-dependent azimuthal correlations relative to spectator 
and participant planes in Au + Au collisions at /Syẹ = 200 
GeV. Phys. Rev. Lett. 128, 092301 (2022). doi: 10.1103/Phys- 
RevLett.128.092301 

[33] X. Zhao, G. Ma, Search for the chiral magnetic effect in col- 
lisions between two isobars with deformed and neutron-rich 
nuclear structures. Phys. Rev. C 106(3), 034909 (2022). doi: 
10.1103/PhysRevC.106.034909 

[34] B. Chen, X. Zhao, G. Ma, On the difference between signal and 
background of the chiral magnetic effect relative to spectator 
and participant planes in isobar collisions at \/5,, = 200 GeV. 
arXiv:2301.12076 

[35] C. Hartnack et al., Modeling the many body dynamics of heavy 

ion collisions: Present status and future perspective. Eur. Phys. 

J. A 1, 151-169 (1998). doi: 10.1007/s 100500050045 

[36] H. Wang et al., Design and tests of the prototype a beam mon- 

itor of the CSR external target experiment. Nucl. Sci. Tech. 

33(3), 36 (2022). doi: 10.1007/s41365-022-01021-1 

[37] W. Huang et al., Laser test of the prototype of CEE time 

projection chamber. Nucl. Sci. Tech., 29(3), 41 (2018). doi: 

10.1007/s41365-018-0382-4 

[38] D. Hu et al., Extensive beam test study of prototype MRPCs 
for the TO detector at the CSR external-target experiment. Eur. 
Phys. J. C 80(3), 282 (2020). doi: 10.1140/epjc/s10052-020- 
7804-2 

[39] X. Wang et al., CEE inner TOF prototype design and pre- 
liminary test results. JINST 17(09), P09023 (2022). doi: 
10.1088/1748-0221/17/09/P09023 

[40] B. Wang et al, The CEE-eTOF wall constructed with 
new sealed MRPC. JINST 15(08), C08022 (2020). doi: 
10.1088/1748-0221/15/08/C08022 

[41] L. Lyu, H. Yi, L. Duan, Z. Zhang, and Z. Xiao, Simulation 
and prototype testing of multi-wire drift chamber arrays for the 
CEE. Nucl. Sci. Tech. 31(1), 11 (2020). doi: 10.1007/s41365- 
019-0716-x 

[42] S. H. Zhu et al., Prototype design of readout electronics for 
Zero Degree Calorimeter in the HIRFL-CSR external-target 
experiment. JINST 16(08), PO8014 (2021). doi: 10.1088/1748- 
0221/16/08/P08014 

[43] Saint-Gobain, BC-408 material. https://www.crystals.saint- 
gobain.com/radiation-detection-scintillators/plastic- 
scintillators/bc400-bc404-bc408-bce412-bc416 

[44] J. Adams et al, The STAR Event Plane Detector. 
Nucl. Instrum. Meth. A 968, 163970 (2020). doi: 
10.1016/j.nima.2020. 163970 

[45] J. Aichelin. Quantum’ molecular dynamics: A Dynamical 
microscopic n-body approach to investigate fragment forma- 
tion and the nuclear equation of state in heavy-ion colli- 
sions. Phys. Rept. 202, 233-360 (1991). doi: 10.1016/0370- 
1573(91)90094-3 

[46] R. Brun et al., GEANT detector description and simulation 
Tool. (1994). doi: 10.17181/CERN.MUHF.DMJ1 

[47] A. M. Poskanzer and S. A. Voloshin. Methods for analyzing 
anisotropic flow in relativistic nuclear collisions. Phys. Rev. C 
58, 1671-1678 (1998). doi: 10.1103/PhysRevC.58.1671 

[48] M. Ding et al., Calibration of the DAMPE plastic scintillator 
detector and its on-orbit performance. Res. Astron. Astrophys. 
19(3), 047 (2019). doi: 10.1088/1674-4527/19/3/47 


[49] J. Adamczewski-Musch et al., Directed, elliptic, and higher 
order flow harmonics of protons, deuterons, and tritons in 
Au + Au collisions at \/8yy = 2.4 GeV. Phys. Rev. Lett. 125, 
262301 (2020). doi: 10.1103/PhysRevLett. 125.262301 

[50] M. S. Abdallah et al., Light nuclei collectivity from ,/syy = 
3 GeV Au+Au collisions at RHIC. Phys. Lett. B 827,136941 
(2022). doi: 10.1016/j.physletb.2022.136941 

[51] P. Russotto et al., Flow probe of symmetry energy in relativis- 
tic heavy-ion reactions. Eur. Phys. J. A, 50, 38 (2014). doi: 
10.1140/epja/i2014-14038-5 

[52] M. Wang et al., Number-of-constituent-quark scaling of ellip- 
tic flow: a quantitative study. Nucl. Sci. Tech. 33, 37 (2022). 
doi: 10.1007/s41365-022-01019-9 

[53] J. Steinheimer et al., Hypernuclei, dibaryon and antinuclei pro- 
duction in high energy heavy-ion collisions: Thermal produc- 


[54] 


[55] 


[56] 


tion versus Coalescence. Phys. Lett. B 714, 85-91 (2012). doi: 
10.1016/j.physletb.2012.06.069 

T.T. Wang and Y. G. Ma, Nucleon-number scalings of 
anisotropic flows and nuclear modification factor for light nu- 
clei in the squeeze-out region. Eur. Phys. J. A 55, 102 (2019). 
doi:10.1140/epja/i2019-12788-0 

T.Z. Yan et al., Scaling of anisotropic flow and momentum- 
space densities for light particles in intermediate energy 
heavy ion collisions. Phys. Lett. B 638, 50-54 (2006). 
doi:10.1016/j.physletb.2006.05.018 

L.M. Fang, Y.G. Ma and S. Zhang, Simulation of col- 
lective flow of protons and deuterons in Au+Au collisions 
at Epeam=1.23A GeV with the isospin-dependent quantum 
molecular dynamics model. Phys. Rev. C 107, 044904 (2023). 
doi:10.1103/PhysRevC. 107.044904 


